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■  SUM-'IARY 

It  is  the  purpose  of  this  paper  to  describe 
an  investigation  of  the  iioriBal  modes  of  lateral 
vibration  of  a  free- free  beam   having  variable  stiff- 
ness and  mass  considering  the  effects  of  bending 
and  shear  deflections  and  rot? ■■■'"•  'nertia,  by  means  of 
an  electronic  system^  A  practical  engineering  pro- 
blem of  this  type  is  found  in  the  vibrations  of  an 
aircraft  v/ing,  or,  as  in  the  case  of  this  investi- 
gation, the  normal  modes  of  vV- '-rlons  of  a  i.aval 

vessel  for  which  data  and  solutions  by  other  methods 
were  readily  accessible „ 

The  components  of  the  electx'onic  system  used 
conslBt   of  standard  feeaii-icic  amplifier  units  whose 
interconnections  are  wlr'^.3.   resistors^  and  capaci- 
tors, relay  switches,  and  associated  pov/er  supplies o 

Results  obtained  considering  bending  deflnrttion 
only  are  in  close  agreement  with  results  obtained 
hj   graphical  methods  and  an  IBM  computer.  Results 
obtained  considering  additional  effects  of  shear 
deflection  and  rotary  inertia  have  IJjuited  accuracy 
due  to  the  assumptions  made  in  setting  up  the  computer 
equation. 

This  investigation  was  conducted  in  the  labora- 
tories of  the  Aeronautical  Engineering  Department, 


tinf-'rersity  of  Michigan ^  the  requirements 

IGT  Master  Degrees,   by  '.h   A,   e  ■'-'--"••■  and  G,   S,   Parrctt 


■  THE  US13  OF  AK  ELECTROT^TG  Al^J^lLOG  COMPUTER 
IN  THE  DETEnillNATiuH  OF  THE  lIoyilvIAL  II0DE3 
OF  LATERAL  VIBRATION  OF  NON-ITHIFOPJ^I  BEMIS 

IKTRODUGTION 

Solutions  of  dynamic  problems  in  fields  otltsr 
than  electrical  have  been  obtained  by  electrical 
analogs  for  some  tiriie^  The  early  use  of  aaiplifiers 
in  the  electrical  analog  netvjorks  added  complica- 
tions due  to  variation  of  tube  characteristics. 
However,  ifith  the  advent  ci  x,^e   stabilized  ieed* 
bsck  eiaplifier  new  techniqu.es  were  developed »  IJith 
these  nevf  techniques  multi-unit  computers  have  been 
designed  in  compact  form  ivhere  the  eleraentJB  are 
irt  erconnected  by  plugging  into  Jacks  at  a  coroinon 
patch  bay.  One  such  comr^iiter  5.3  the  Rneves  Elsct,:-onic 
Analog  Computer  developed  by  the  Ueevea  Instrument 
Corporation  for  the  Special  Devices  Center  of  th® 
Office  of  Haval  Researchf 

The  techniques  employed  herein  are  based  on 
individual  stabilized  feedback  amplifier  units 
.designed  after  the  circuits  given  in  Ref ,  4?  hy 
C.  E.  Hox-^e  and  R.  Howe,  for  their  invest igatio:.is  vrith 
electronic  computers  conduated  at  the  University  of 


6 

Michigan  during  the  sur  The  : 

Investigations  has  not  yet  been  pabli.'^hnd.  Th^:  m\i±n^ 

ment  developed  during  the  above  invsstigatl: 

;32ade  available  to  the  authors  for  their  ur^^   -in  this 

project « 

This  equipment  cor.si --''^^  r,r   r,?vr-^-li^  fr-rib^ck 

amplifier  units  which  by  simple  external  charges 

in  connections  v^ill  serine  as  i.,.....,,-....:ooi:',  '- r- 

entiator^  multiplier  or  i:  slgii 

Because  of  this  ability  to  perf cm  r.ho  op     one. 
of  aritlimetic  and  calculus  on  an  input  volt     /.e 
units  are  temed  "operational  amplifiers."  Provided 
the  associated  power  s     33  a.fe  well  re         d 
the  circuit  components  are  of  higli  precision ,  the 
a?c\iracy  of  the  results  :  at  ions  is 

adeuqate  for  niost  engine:  ^'^  *  problems,  Th®  com- 
ponents of  the  computer  other     the  operational 
aiaplifiers  and  their  associ-.-^v.  ^..,..^:r  supplies  con- 
sist of  'fd.r&B,   resistors,  capacitors,  and   '  ,y 
svritches  for  use  in  ssuo  oonditiors  and 

introducing  tlie  variables  of  '  ^' 

systea.  Solutions  in  the  xc^.m  cf  output  voltagee 
are  amplified  and  recc.    on  a     r  tape  rtcoi 

The  Hov-res*  investie^^^^^i'i  i     ^   solutions  of 
problems  dealing  with  de       is  and  norEr.^l  modes 
of  vibration  of  beams  for  various  loading  and  end 


conditions,  A  solutioxi  of  the  natural  modes  of 
vibration  of  a  free- free  beam  v/ith  variable  moments 
of  inertia  and  variable  mass  was  obtained  con* 
sidering  bending  deflections  onlyc  The  example 
used  V7as  the  hull  of  a  ficticious  ship  whose  moment 
of  inertia  and  mass  distributions  vrere  symmetrical 
about  its  middle  section. 

The  primary  problem  concerned  herein  ia  a<n 
extension  of  the  above  in  that  moments  of  xx.icr'&ia  and 
mass  distributions  are  not  izymmetrical  aboi't  the 
middle  section  and  the  ef''"  '  ■■?.  of  shear  dei.;.ccticnB 
and  rotary  inertia  are  ihcluded.  Data  for  the  problem 
imre   taken  firom  an  actual  Tessel^  the  APA  Hy^   of  the 
United  States  Navy, 

An  investigation  of  the  noiiiial  modes  of  vibra** 
tion  of  the  same  vessel  by  anal;<7t:lcal  method!?.  (Stadola) 
has  been   conducted  by  the  University  of  Michigan, 

Engineering"  ITechanics  Deparlyment^  in  conjiinction  vdth 

•7 
work  under  contract  for  the  United  States  Navy,  Solu- 
tions were  obtained  for  the  first  and  second  nomal 

t. 

modes  considering  bending  deflections  only^  Under 
another  United  States  Navy  contract  to  the  International 

tesB  Machines  Corporation,  solutions  of  the  first 
five  normal  modes  considering  bending  deflect icns  only, 
and  also  id.th  additional  effects  of  shear  deflections 
and  rot£.ry  inertia  were  obtained  by  means  of  a  mechanical.- 
electrical  computing  machine*  Results  obtained  by  the 


above  tv70  methods  are  cQiaparad  with  those  obtained  in 
the  follovdng  investigation  wJilch  uses  a  relatively 
simple  electronic  analog  coapu-oer. 

Several  simpler  practicdl  problems  dealing  with 
aeaas  are  first  solved  in  oraer   to  familiarise  the 
authors  v^ith  the  use  of  the  equipment  and  to  demon-' 
strate  the  simplicity  of  computer  set-ups  ior  5.n« 
cj'easing  complicity  of  the  differential  equaticiis 
to  be  solved „  These  problems  are  followed  b]*'  the 
solutions  of  the  first  three  narm^l  modes  of  vibration 
of  the  APA  &'7   considering  bending  deflections  onlyj. 
and  of  the  first  four  normal  mode;.-:  considering  the 
additional  effects  of  shear  deflections  and  3*otary 
inertia , 

The  authors  •vd.sh  to  express  ^neir  appreciation 
to  Associate  Professors  E,  J.  Les'ier^  M,  H,  WicholSj 
and  L,  L,  Rauch^  Aeronautical  Engineering  Department, 
University  of  Michigan,  for  th-?'1,r  assistance j,  coopera-- 
tion,  and  guidance  in  carrying  out  this  investigation. 


EQUIPiMMT  AND  FROCFDUF^ 

As  stated  in  the  introduction  the  various 
el  Clients  necessary  for  setting  np  a   oosnplete  computer 
circuit  were  previously  assernblec.  Exid   available  as 
laboratory  equipments  A  complete  description  of  the 
elements  are  to  be  found  in  Rex .     oince  tnis 
reference  is  at  present  t5.me  rinpublishedj  brief  de- 
scriptions of  the  various  components  used  are  given 
herein^ 

The  unit  upon  which  the  computer  is  based  is  a 
direct  current,  thres-stn.q-f?^  feedback  amplifier 
having  good  stability,  high  gain  and  an  effective 
phase  shift  of  ISO  degrees.  In  the  following  pro- 
blems the  amplifier  is  uesd  as  a  multipler  or  divider ;, 
sign  changer 5>  and  integrator.  To  have  the  amplifier 
perfoim  any  of  these  operations  it  is  necessary  only 
to  make  external  changes  of  the  impedances  in  the 
input  and  feedback  circuits.  Fig.,  1  shows  the 
circuit  diagram  of  the  operational  amplifier-  Fig«  2 
is  a  photograph  of  the  chassis  in  which  the  circuit 
is  mounted.  Fig.  3  and  Fig„  4  are  photogx'aphs  of 
the  operational  amplifier  as  r.njltiplier  and  integrator., 
respectively.  Due  to  the  11^0   deg-voes  phase  shift  in 

* 

the  amplifier  it  changes  the  sign  of  the  input  volt- 
age in  every  operation*  For  pure  sign  changing  the 


input  and  feedback  impedances  consist  of  equal 
resistors;  for  multiplication  the  desired  ratios 
of  feedback  to  input  resistors  are  plugged  into 
the  amplifier  circuits 

The  amplifiers  are  connected  to  a  power  supply 
distribution  panel  by  means  of  st;:-v7irG  shielded 
cables.  The  txiro  knobs  shorn  on  the  chassis  of  the 
amplifier  are  used  to  balance  the  amplifier  for  zero 
direct  current  output  prior  to  usi  in  the  computer 
circuit.  The  knobs  are  connected  to  two  variable 
resistors  associated  vdth  the  inptit  tube<, 

To  balance  a  multiplying  or  .>ig*n  changing  ampli- 
fier, the  input  is  shorted  to  ground  and  equal  resistors 
(one  megohm)  ere   plugged  into  the  input  and  feedback 
circuits.  By  use  of  a  multi-rang?  direct  current 
vacuum  tube  voltmeter  connected  butxireen  ground  and 
the  output,  zero  output  is  obtainod  by  adjusting  the 
two  knobs. 

The  procedure  for  balancing  the  integrating 
amplifier  is  similar.  Uith  a  one  megohm  resistor 
and  a  one  oocrofarad  capacitor  as  input  and  feedback 
inlpedances  respect ively,  balance  is  obtained  by  adjusting 
the  knobs  for  a  constant  output. 

The  procedure  for  balancing  combinations  of 
amplifiers  in  a  computei-  set-up  will  differ  sosaewhat 
depending  on  the  problem  or  combirationo  A  small 


unbalance  in  an  individual  amplifier  may  build  up, 
particularly  through  integrating  amplifiers,  into 
unacceptable  overfall  unbalance  of  the  computer,  A 
procedure  for  balancing  a  complete  computer  circuit 
of  six  amplifiers  is  described  later  in  the  problem 
dealing  v;ith  the  vibrations  of  a  ship's  hull. 

The  jacks  seen  in  the  chassis  of  the  amplifier 
were  designed  for  banana  plugs  and  spaced  for  General 
Radio  Type  274-H  double  plugs.  Input  and  feedback 
resistors  are  Continental  X~type,  accurate  to  p3.uf3  or 
minus  one  percent.  Polystyrene  capacitors  used  in 
the  feedback  circuits  of  integrating  amplifiers  ar© 
Western  Electric „  one  mocrofar.id  capacitors  having 
high  leakage  resistance  anu  low  dielectric  absorption o 

The  high  voltage  povfer  supplies  required  for  the 
aiflplifiers  consist  of  plus  and  minus  350  aiid  minus 
190  well  regulated  and  filtered  direct  current  volta.'^es. 
Fig,  1,  These  voltages  together  '.-dth  a  six.  volt 
supply  for  the  amplifier  heaters,  are  taken  to  a 
distribution  box  by  means  of  shielded  cables «  The 
distribution  box  used  has  twelve  outlets  for  the  con--- 
nsctions  of  amplifiers.  Although  both  alt cths ting 
current  voltage  and  direct  current,  voltage  eripnlies 
were  available  at  the  distilbutioii  box  for  the  ampli- 
fier heater  circuit j  the  latter  v^as  found  to  reduce 
the  amount  of  sixty  cycle  oscillations  which  wer© 
Gxpeoially  noticeable  in  the  computer  outputs  when 


several  integrating  ?5mpli,fisrs  trerein  use. 

Due  to  the  low  output  impedance  of  ths  operational 
sxaplifiers  the  output  current^  is  small  necersltrting 
a  direct  current  power  amplifier  with  a  high  ratio  of 
input  to  output  impedances  in  order  to  obtain  satis*« 
factory  records  of .  the  computer  outputs*  For  this 
-.-^rpose  Brush,  Model  BL-.913,  direct  current  amplifiers 
were  placed  between  the  desired  outputs  o-:  the  computer 
and  the  recording  oscillograph*  "lUhe  TQcardiim   op.eillo- 
graph  used  i/ms  a  two  channel,  Bruuh,.  Model  jj^y^dOZ^ 
■ .   :netic  typ®* 

The  equipment  used  for  simulating  the  variable 
^'"-'efficients  of  the  equatione  and  setting  the  ini- 
tial conditions  of  the  problems  consists  of  a  syn- 
chronous contactor,  un-^,-^?  of  stepping  rel^^'-  "^^^tehes 
%lth  their  panels  for  lug- in  resistors,  initial  condition 
relay  switches,  and  a  s^elay  control  panel* 

Fig.  5  shows  a  photograph  of  a  stepping  relay 
unit  with  panel  having  the  proper  lug- in  resistor 
assemblies  for  simulating  the  variable  moaests  'of 
iixcii-t<ia  along  the  length  ..-.  a  shipVs  hull  used  in  a 
later  discussion.  The  stepping  relays  are  dri¥es^ 
^^   ei(if;trie  pulses  from  a  synchronous  mo'     .-  .^w^i^i.^ 

synchronc         Iriv©      ft  at 
ons  re^9-0i.u-€>j.on  per  secor  ' ..  ju   vu-y  oi:.cv.LO     --^'unted 
two  cams,  one  having  one  flat,  th-s  other  having  four 


flats*  An  tlnde  Lab.,  Inc.,  typo  32-2RLS  microswitch 
rides  on  each  of  the  cams;  the  one  giving  four 
pulses  per   second  drives  the  stepping  rela-ys;  the 
si^itch  giving  one  pulse  per  second  is  sent  directly 
to  the  relay  control  pansl  for  th-3  purpose  of  start- 
ing the  computer  probl^^nis  alvjays  on  the  same  flat  of 
the  stepping  relay  cm »     This  unifonnity  is  provided 
to  minimis©  the  effects  of  the  machining  tolerances  of 
the  cam  flats  on  the  timing  of  the  probleni« 

The  stepping  relay  circuit  diagram  is  shown  in 
Fig,  6.  Each  relay  consists  of  three  arcs  of  forty 
contacts  each.  Two  of  the  contact  ares  have  bridging  . 
wipers  while  the  third  has  a  non-c ridging  'tapei". 
The  latter  is  used  to  impose  and  iremove  initial  con- 
ditions V7hile  one  of  the  former  is  used: to  cv;t  in 
the  series  connected  resistors  on  the  lug«in  panel* 
With  this  equipment  it  is  therefore  possible  to 
simulate  a  variable  function,  of  a  beam  for  instance „ 
in  forty  stations  along  the  beam.  The  number  of 
stations  used  together  vjith  the  number  of  relay  contacts 
wiped   per  second  determine  the  length  of  the  problem 

beam  in  terms  of  seconds.  To  obtain  a  correct  solutiotj 
to  a  problem  then,  the  imposed  end  conditions  must 
be  satisfied  vjithin  the  exact  length  of  the  beam  as 
determined  by  the  stepping  relay  panel  set-ri:  and 
motor-driven  cam  combination.  Such  a  combination  is 


described  later  for  the  varying  properties  of  a 
ship^s  hull  which  are  varied  in  twenty  discrv^tc  steps 
along  the  length, 

A  stepping  relay  control  panel  circuit  designed 
to  control  three  such  stepping  relays  plus  the  imposing 
and  removing  of  the  initial  conditions  is  shcxm  in 
Fig,  7*  The  follovjing  description  of  thivS  circuit  is 
talc  en  directly  from  Ref «  6, 

Relay  F  is  the  master  pulsing  relay,  its  pulsing 
rate  depending  upon  the  two  cains  on   the  s'ynchronous 
contactor  and  on  the  position  of  a  remote  &>tdtcha 

Relay  G,  through  normally  closed  contacts,  passes 
pulses  from  relay  F  to  the  coil  of  stepping  relay  A« 
1/hen  stepping  relay  A  reaches  position  40,  relay  G  is 
energized  and  no  longer  passes  pulses,.  Stepping  relay 
A  then  stops.  Relays  H  and  J  perform  the  ssme  functions 
for  stepping  relays  B  and  C, 

These  three  relays  G^  H  and  J  also  play  an 
important  part  in  5.mposing  the  initial  conditions, 
^'Jhen  all  three  of  these  relays  are  energised  {vfnen 
all  three  stepping  relays  are  on  contact  40)  power 
is  furnished  to  the  coil  of  relay  Lj  which  is  then  closed, 
This  remo'ves  power  from  the  "locking^'  contacts  on  relay  M, 

Relays  L,  M  and  H  perfonn  the  functions  of  auto- 
matically imposing  and  removing  the  initial  conditions » 
The  initial  conditions  are  imposed  as  soon  as  all  three 


Xi. 


stepping  relays  reach  contact  40, 

•/Then  relay  N  is  closed ^  the  initial- condition 
relays  are  energized j  thereby  remc^/lng  the  initial 
conditions.  Relay  N  is  controlled  by  normally-open 
contacts  on  relay  M«  If  relay  M  is  momentarily  ener- 
gized it  remains  closed  by  irirtue  of  its  "electrically 
locking"  contacts.  These  contacts  obtain  power  from 
normally- closed  contacts  on  relay  L,   (As  long  as 
relay  M  is  closed,  relay  K  is  closed  and  all  initial 
conditions  are  removed.)  If  relay  L  is  energised  (all 
stepping  relays  on  contact  40)  relay  M  "drops  out" 
and  the  initial  conditions  are  restored.  The  initial 
conditions  are  not  reraoved  raitil  relay  M  is  again 
energi'^ed  which  is  done  as  scon  as  any  one  of  the  stepDlng 
relays  readies  contact  1„ 

The  stepping  relays  always  stop  on  contact  /.-.n.. 
t/hen  they  are  in  this  position  relays  G,  H  and  «;  are 
energized  and  no  longer  furnish  drixdng  pulses  to  their 
respectiTe  stepping  rela^rSo  Relay  L  is  ensrgizedj 
renao-ying  poii/er  from  the  '^^loc.<:ing"  contacts  of  relay  M. 
Relays  M  and  N  are  inoperative,  no  power  is  furnished 
GD  the  initial- condition  relays  and  the  initial  conditions 
are  imposed. 

Relay  0  is  the  starting  relay »  controlled  by  the 
remote- control  momentary- contact  starting  butt on g  S^  • 
XJhen  this  s;^itch  is  closed  momentarily,  relay  0  is 
ener'^ised  as  soon  as  the  nex-t  pulse  is  furnished  by  the 


microswitch  on  the  one-per« second  cam.  Relay  0  will 
then  renain  closed  until  rela3/'  M  closes «  As  soon  as 
contact  40  is  left,  relays  F,  G  and  H  open  and  pulses 
are  continued  to  be  supplied  to  the  stepping  relays. 
Simultaneously  relay  L  drops  out,  energising  the 
"locking"  contacts  of  relay  M«  At  the  instant  any 
one  on  the  stepping  relays  reachers  contact  1  relay  M 

closes  and  remains  closed.  This  immediately  remoYes 
the  initial  conditions  and  deenergises  the  starting 
relay  0. 

In  ease  any  one  or  more  of  the  three  stepping 
relays   are  not  used  the  corresponding  switches  S^j, 
Sq  or  Sc  should  be  closed*  This  will  then  permit  noiiaal 
operation  of  relay  L* 

Figo  S  is  a  photograph  of  a  jcmplete  computer 
set  up  to  solve  a  fourth  order  differential  oquation 
with  variable  coefficients.  Fig.  9  is  a  box-type 
diagram  of  the  same  set-up  identii'ying  the  -various 
components  of  the  systsra»  Diagra;;.i3  shov/ing  one  arc-^^ 
temal  circuits  of  the  amplifier  3ombination3  for 
performing  the  particular  operations  reqiiirea  are 
given  in  the  foilos^ling  discussion  of  the'  iridivic 
proble7iis. 
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DISCUSSION  AMD  RESULTS 

Preliminary  to  the  primary  object  of  this 
investigation,  there  is  prssented  the  solutions 
to  tv7o  simple  beam  problems  as  determined  hj  the 
analog  computer.  Such  detemiinatioas  were  made  as 
part  of  a  program  of  familiarization  with  the 
technique  to  be  used  and  to  develop  a  facility  in 
setting  up  and  operating  the  ecmponents  ox  Gi.e 
electrical  circuits  which  become  integral  parts  of 
the  computer  netvjork  used  in  the  solution  of  the  main 
problem. 

Pai't  j: 

The  first  preliminary  problem  was  the  dGtemin*- 
ation  of  the  static  deflection  under  uniform  load 
of  a  horizontally  supported  beam  of  constant  crosfj 
section  ^-fhich  is  small  in  comparison  with  its  len,^h» 
Two  different  types  of  end  fixity  are  demonstrated; 
■^~is .  5  clamped  and  hinged « 

The idiff ©rent ial  equation  of  the  elastic  curve 

9 

of  such  a  beam  is  given  as 

FT   2JL  r   vv()c)  .  (1) 


where:  w(x)  is  the  load  per  unit  length  along  the 
beam* 
y  '•  vertical  deflection  of  the  bG^m  at 

any  point  x, 
X  «  distance  e.long  the  beam  measured  from 

one  end, 
E  =  Young's  Hodulus  of  Elasticity « 
I  «»  area  moment  of  inertia' of  a  C3"oss 
section  of  the  beam  'dth  respect  to 
the  centroidal  axiSa 
Bending  moment  and  shear  force  at  any  point  x 
along  the  beam  are 


Banding  Moment :  M  /^x  --  p  \ 


Shear : 


Q^x^  --t 


X. 


d  X- 


Type  I 


Beam  Clamped  at  Both  Ends 


\  Efltf" 


End  Conditions:  Zero  slope  and  aero  deflection  at 

each  end.  These  boundary  conditions 

are  expressed  as 

y(0)  =  yHO)   '^   y(L)  »  yML}  '^  0. 


-^J' 


9 
Theoretical  Solution: 

y(x)    =    ilM  i±  ^   ..  fa3   ...  x'^  ^ 

12  EI  V2  2    ; 

y(max)  =  iiLML'      @  X  «  L/2o 
364  EI 

The  computer  equation  is  set  up  by  making  a 
change  of  the  independent  variable  in  the  original 
equation  (1}  ..  The  independent  variable  x  is-  changed 
to  t,  time  in  seconds,  and  the  length  of  the  original 
beam  is  expressed  at  Tj  total  elapsed  time  for  solu«» 
tion  in  seconds » 

Thus  -  -  L  t  J  and    d""  {  )  ^„  f    d"  (   ) 


dx*^     T''  dt"^ 

The  ebmputsr  equation  is  then 

EI  d^I  „  L^  -^  {^) 
dt^   T^ 

(2) 

The  computer  circuit  for  the  solution  of  this 
equation  is  given  in  Fig.  10,  The  end  conditions 
j{Q)   =  y*{C)  ==0  are  obtained  by  initially  shorting 
the  feedback  capacitors  of  A^  and  A,^  Shear  force 
and  bending  moment  have  definite  valuer  at  the  ends 
of  the  beam  and  are  simulated  by  battery  voltages 
-V'a  and  Yb   respectively  initially  applied  to  the 
(n|"      capacitors  of  Ai  and  Ao.  As  the  va3.ues  of  shear 
!J^      force  and  bending  moment  are  unknuwiSj,  so  the  applied 
voltages  -Va  and  Vb  to  be  applied  are  at  first  unknoi-Jii. 


Vb  vms  made  fixed  at  about  sir.  volts  and  -Va  was 
varied  until  a  correct  solution  was  obtained,  A 
constant  battery  vo3.tage  V  was  applied  as  input 
voltage  to  Ax  to  simulate  the  uniform  loading^  w(x) ^ 
of  the  beam  J  and  was  mea  spired  in  terms  of  recorder 
deflection  units.  The  outputs  of  A/^  and  A3  were 
connected  through  amplifiers  to  channels  1  and  2  of 
the  Brush  recorder  for  recording  oscillographs  of 
-y'  and  y« 

A  correct  solution  was  obtained  in  the  folloifing 
manner:  With  a  constant  i23.put  voltage  V  of  about 
13   volts  applied  to  A^j,  Vb  constant  at  6  volts,  and 
-Va  set  at  an  arbitrary  value ,  all  initial  conditions 
were  imposed  by  closing  the  initial  condition  relay 
switches.  The  problem  was  started  by  de- energizing  all 
the  initial  condition  relays  simultaneously  thus 
releasing  the  end  conditions*  Several  trials  were 
made  with  different  settings  of  the  potentioineter 
controlling  the  battery  voltage  '•Va  before  the  proper 
end  conditions  ^irere  satisfied  and  recorded  on  the 
oscillograph.  Fig.  11  shovjs  the  solution  of  the 
problem  as  recorded.  The  length  of  the  solution,  T, 
on  the  oscillograph  is  T  =  3.6  seconds.  V  »=  1«3.9 


i.7 


Then  wM  L^  ==  VT^  -  1»3  x  l67o96  =  21^.2 


y(max)  theoretical  ==  2_lg_>2  -«  0»56&\ 
yfmax)  from  Fig.  11       «  Oo55Q. 

T^ge„I_X     Beam  Hinged  at  Both  Ends 


End  Conditions:  Zero  deflection  and  sero  bending  monient 

at  each  end.  These  boundary  conditions 
are  expressed  as 

y(0)  =  y^MO)  =*  y(L)  «  y«»(I')  -^  0« 

g 

Theoretical  solution:  ^ 

y(niax)  ==  _5   ^'/(2^_)  L^      @  sc  -  L/2. 

q^/     EI 
The  computer  equation' xs  the  same  as  fo2»  Type  I, 

but  a  change  is  necessary  in  the  computer  circuit 

due   to  the  change  in  end  conditions.  The  coraputer 

circuit  for  the  solution  of  this  problem  is  gi"/en 

in  Fig.  12.  The  end  conditions  yfO)  »  y»«(0)  «  0 

V7ere  obtained  by  shorting  the  feedback  capacitors 

of  A2  and  ki^o     Shear  force  and  slope  have  definite 

values  at  the  ends  of  the  beair:  and  are  simulated  by 

battery  voltages  -Va  and  Vb  respectively  applied  to 


capacitors  of  Ai  and  A3,  The  correct  solution  is  ob- 
tained as  before  for  Type  I.  Several  trials  again 
were  necessary  in  varying  the  potentiometer  of  -Va 
to  obtain  the  proper  end  conditions. 

Fig,  13  shows  the  solution  of  the  problem  as 
recorded  on  the  oscillographs  ^^he  length  of  the 
beam  l\   on  the  oscillograph  is  measured  at.  T  "  2.% 
.3eco;nds.  V  «  1,3 » 

Then  w{x) lA  «  VT^  -  1.3  :?£  76o?  «  99.6 

EI 

y{ma:x)  theoretical  ~  5  -^  99.6  »  1,30 
y(iaax)  from  Fig,  13         »»  1.15 

The  second  preliminary  problem  was  the  deter- 
mination of  the  first  three  noimal  modes  of  lateral 
vibration  of  a  uniforai  free- free  beam  considering 
the  effects  of  bending  deflection  only.  The  vibrating 
heem.   is  considered  loaded  by  inertia  forces  due  to 
its  Q\m  mass  and  acceleration. 

The  differential  equation  of  motion  of  the 

;n 
elastic  curve  of  such  a  he^a   is  given  kxj   ■' ' 


where  /^  «  A  y  3  the  mass  distrib^ition  aloruj  ohe  beam 
f  «  the  density  of  the  material  of  the  beam. 


A  =  the  cross  sectional  area  of  the  beania 
g  «  the  acceleration  due  to  gravity* 
//'^('^•^^   the  inertia  forces  acting  on  the  beam« 

y  *  vertical  deflection  of  the  beam  at  x» 
X  =  distance  along  the  beam  measured  from 

one  end, 
E  =  Young's  Ilodulus  of  Elasticity „ 
I  «  Area  moment  of  inertia  of  a  crross 
section  of  the  beam  x-rith  respect  to 
the  centroidal  axis^ 
It  is  assumed  that  yCrs^t)  =  X'Cx)  e*^"""*"^ 
VJliere  X(x)  is  a  function  only  of  Xj  and  ir,  independent 
of  timosf  e^'*  ^  represents  sinusoidal  oscillations  of 
frequency  CO. 

Then      .]£y{3c,t)  =  =^X(x)u)^s^*"^^ 
and  equation  (3)  becomes 

dx^ 

or     ;, 

EI  d'^'X  -  X  -  0 

iU) 
The  computer  equation  is  set  up  as  before  by 
making  a  change  of  the  independent  variable  in  the 
original  equation  (4) .  The  independent  variable^  x^ 
is  changed  to  t,  time  in  seconds,  and  the  length  of 
the  beam  is  expressed  as  T,  total  elapsed  time  for 
solution  in  seconds* 
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Then  x  =  L  t  and   (j  (  )  •=  L^  cf{   ) 
The  computer  equation  becomes 
EI  T^  d^X  -  X  «=  0 

For  simplicity  we  let 

2 


C5) 


EI 


and 


the  natural 


frequency  of  vibration  for  the  n"^^  raodeo  In 
addition,  for  the  computer  equation  (5)  we  denote 

so  that  the  computer  equation  reduces  to 


C  d% 
dt^ 


X 


(6) 


For  the  computer  solution g  the  C  was  given  a  value 
of  unity  (1  megohm)  and  the  problem  was  solved  by 
firming  a  lengthy  Tj  on  the  oscillograph  solution  for 
which  the  simulated  end  conditions  as  determined  by 
the  beam  supports  were  meto 


-^^ 


.-L 


The  computer  circuit  for  the  vsolution  of  equation 
(6)   is  giiren  in  Fig.   lU,     The  end  conditionv'3  to  be 
satisfied  for  a  free-free  besja  are  that  the  bending 
nioraent  anij.  shear  forco  at  each  end  are  zero.,     These 
boundary  .'^oMitions  are  expressed  as 

X"()l    =  X«««(0)   «  X«»(L)    «  X««ML)   «  0 

Tj>  satisfy  these  end  conditions  on  the  computerj, 
the  .feedb.K-k-  capacitors  of  A^  and  A2  are  initially 
shor.'teda  /s  there  is  a  definite  but  unknoKH  slope 
and  deflejJ  ion  at  each  end  of  the  beanij  these  are 
sir/.ulated  on  the  coiaputer  by  battery  voltages  -Va 
arjl  ?b  rer'pectively  initially  applied  to  the  capaci- 
to7.'3oiA^   ar*i  A/j.o  As  beforf^,  Vb  was  fisced  at  about 
sIj:  YOlts  and  «Va  i7as  varied  for  different  trial 
so*i,-;tion£i/  ^mtil  the  end  conditions  of  se:"o  shear 
f or  ;e  and  ioending  moment  were  satisfied ^ 

The  jut puts  of  Ai  and  A2  ware  connected  through 

sHii'lifier**  to  the  two  channels  of  the  Brush  recorder 

ior  recopjing  oscillographs  of  X'^  and  -X'^^o 

jorrect  f;lutions  showing  the  fulfillment  of  the  end 

conditicfi I  required  v/ere  obtained  v^^hen  the  minimum 

or  maxir  jir of  X^\   depending  upon  the  number  of  the 

mode  J  pc.Miied  through  the  zero  axis.  The  -X"*  curve 

was  used/  In  measuring  the  lengthy  Tp  of  the  solution. 

I 
This  cu;'(  e  was  used  in  preference  to  X'^  because  the 

^X»^'  cri'vc  has  a  definite  finite  slope  at  each  end 

of  the  stjlution. 


Correct  solutions  of  the  frequency  for  the  first 
tm^ee  nonnal  modes  of  vibration  of  the  beam  were  ob- 
tained in  a  manner  similar  to  that  previously  des- 
cribed. Vb  was  made  constant  and  -^Va  was  initially 
set  at  an  arbitrary  valua*  All  initial  end  conditions 
were  imposed  by  closing  the  initial  condition  relay 
switches.  The  problem  was  started  by  simultaneously 
releasing  all  the  end  conditions.  Sever?!  trial 
settings  of  the  potentiometer  controlling  the 
Toltage  -Va  were  necessary  before  a  coi*i'ect  solution 
for  each  mode  v/as  obtained. 

Solutions  of  the  first  mode  V7ere  quite  readily 
obtained,  but  for  the  second  and  third  modes  the  set-» 
ting  of  the  potentiometer  controlling  the  -j/oltage  -Va 
was  found  to  be  very  criticalo  The  inherent  slight 
instability  of  the  amplifiers  used  and  v'nriations  in 
power  supply  voltage  V7ere  enough  to  cause  trouble  in 
repetition  of  solutions.  I-feny  trials  tvf.-ars  necessary 
to  obtain  a  few  correct  solutions.  Fig,  15  shows 
the  correct  solution  oscillographs  obtained  for  tfce 
first  three  modes.  The  results  obtained  for  a^^  agg 

btA  aj   checked  very  closely  vdth  those  given  in  thf5 
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Ajjpendix  of  Den  Hartog  '  i 


Mode 

Den  Hart^ 

Coraguter 

Otj 

1 

22,4 

22e4 

aa 

2 

61,7 

•• 

Gim 

«3 

3 

121,0 

121,0 

:^3 


In  ©Taiuating  the  results  of  the  preliminary 
problems,  several  factors  were  noted  as  being  most 
important  in  the  obtaining  of  correct  solutions 
on  the  analog  coinputer  of  the  laboratory  type  used. 
These  factors  are  mentioned  briefly  here  znd  were 
kept  constantly  in  mind  in  setting  up  and  operating 
the  ccmputer  for  the  main  problem. 

1*  The  power  supply  to  the  individual  amplirit^rs 
must  be  as  required  by  the  amplifiers  and  as  non- 
variant  as  practically  possible.. 

2^  The  tube  heater  voltage  of  about  6  voltjj 
is  best  obtained  frora  a  stOrag-^  battery. 

3.  Each  indi'^ridual  aj:]iplifisr  should  be  finely 
balanced  for  zero  gain  and  this  balance  .should  be  kep'i 
by  frequent  check, 

4.  The  amplifiers  should  be  finely  balanced  for 
sero  gain  in  the  eompuber  network  in  groups  of  2,  3* 
and  4,  and  be  kept  in  balance  by  frequent  check, 

5,  All  precautions  should  be  taken  that  the 

60  cycle  '^pick-up"  by  the  notvrork  be  kept  to  a  minimum 
that  can  be  tolerated  by  the  system, 

6,  As  the  measurement  of  length  of  time  as 
recorded  on  the  oscillograph  is  Diost  critical  in  arriving 
at  a  solution,  pen  lag  of  th"e  recorder  should  be 

kept  at  a  minimum,  and  the  110  volt  A,C.  60  cycle 


driving  the  syuchronous  motor  of  the  BriAsh  recorder 
should  be  carefully  regulated.  Pen  lag  effect  can 
be  effectively  reduced  by  increasing  the  ariT.plifier 
output  voltage  recorded  on  the  cscillograph» 

7<.  The  potentiometer  controlling  the  battery 
voltage  «?a  had  three  degrees',  of  fineness  of  control 
which  vfas  found  quite  necessary  for  solutions  of  the 
higher  modes.  Va  for  the  ejeeond  mode  \mn   found  to  be 
very  close  to  the  Va  for  the  third  modeo 

It  was  felt  that  the  accuracy  of  th^j  results 
obtained  in  the  preliminary  ifork  was  well  within  the 
accuracy  of  the  laboratory  t3rpe  equipment  usede  As 
the  technique  on  the  part  of  the  operators  of  the  equip--^ 
ment  improved  both  in  setting  up  the  problem  and  in 

checking  the  balance  of  the  amplifisrs  in  the  netVTorkj 
the  results  obtained  by  the  computer  definitely  improved 
in  accuracy. 

Part  II 

For  the  purpose  of  this  study >  the  vibrating 
ship  is  replaced  by  an  ideal  floating  beara  with.  uni-= 
form  clastic  and  ineitia  properties.  The  differential 
equation  for  the  vertical  motion  of  the  elastic  curve 
of  this  floating  beam  has  been  derived   including 
the  effects  of  bending  and  shear  deflection ^  rota:ry 
inertia,  external  loading ^  damping  force ^  and  buoyancy 


force.,  The  problem  here  considered  is  the  solution 
of  the  differential  equation  for  the  vertical  motion 
of  the  elastic  curve  of  the  vibrating  ship  considering 
first  bending  deflection  only^  and  ser.o--r^,.  'hf^rding 
and  shear  deflection  *v^'ith  rotary  inertia  effect* 

The  complete  differential  equation  of  the  elastic 

.  11 
curvQ   IS  given  by 


^ 


t2 


r        2 


'') 


Y'   EI  A 


0 


(7) 


where  each  term  has  a  physical  interpretation 

The  first  term 

2  /  .,      .2 


is  the  sum  of  the  moments  due  to  elastic  defosr^ation 

in  bending  and  rotarjr  inertia.     The  differential 

V  2 

opejrator    j.        reduces  th.ls  to  an  equivalent  diistri- 


buted  load, 

The  term  //  H, JLfe '^  ^   -^^  '^'^■^®  distributed  load  du© 

t 


o-f" 


to  translatory  inertia. 

Considering  first  the  effect  of  bending  deflection 
only,  the  differential  equation  reduces  to 


K/  '%J 


1 


2 


{^) 


ifhere  //=  the  mass  distribution  along  the  beam  and 
includes  the  virtual  mass  which  is  an 
equivalent  mass  added  to  that  of  the  ship 

r 

to  represent  the  inertia  effect  of  the 
water  accelerated  \irith  the  snipes  vibra- 
tion. 
For  sinusoidal  oscillations  of  frequency  a>  ,  it 
is  assumed  that  jix^t)   ~  X{x)   e*^"^ 

where  X{x}   is  a  function  only  of  x  distance  along  the 
beam  and  is  independent  of  time„ 
Qjfl^t  represents  sinusoidal  oscillations  of 
frequency  od  <, 


X(.x) 


o  A- 2 
and  assuming  that  E  is  a  constant  and  that  I  andyU 


are  functions  only  of  x  along  the  beam,  equation 

fS)  bscomes 


,2 


dyj 


EI  d2x 


2 


(9) 


vie   let  I  "^^  IgiCx)  and^  ==••  ^q^(x)  inhere  I^  and 
^/ Q  STB  maximum  values  of  moment  of  inertia  and  mass 
respectively „ 

To  set  up  the  computer  equation j  a  change  in  the 
independent  variable  is  again  necessary «  IJhere  x 
i\ras  the  independent  i^ariableg  let  t,  t:^-me  in  seconds 
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on  the  computer  solution  be  the  neu  variable,,  Then 

t  "  ?  K 

L 

where  L  is  the  length  of  the  ship  actually,  T  is  the 
length  of  the  ship  on  the  computer  solution,  !Jhere 
0  fs  X  <^   L  before,  pow  0  &  t  ^  T  for  the  computer « 
In  general^    j'^,  .  .^^^     ^-ti,  . 

and  equation  (9)  becomes 


T^-  d^  /eI  d^X  )  ,, 

?^d^      -^r/' 


(jO 


'X   ==  0. 


(10) 


Substituting  for  I  and^  j  E  «  a  constant,  and 

2 
dividing  through  by  a)     j   and//    ,  we  have 


T^  EI|,     d^ 
yCfo  lV'  df^  ^ 


i(t)   d% 


dt'- 


p{t)X  »  0 


(11) 


Letting o(,^  °  >yo^i  L^         and  C  -  T^ 


EI 


o 


7? 


o( 


'/V\ 


the  natural  frequency  of  vibration  of  the  '!*'•'*  mode 


Oi 


<eii»,  iiii»  w«  nm— 


The  computer  eolation  becomes 


C  d' 


dt^ 


i(t)   d%  p(t)   X  =  0 

Ml  Mill  I  I  I  I  «X'  «  ^  * 


(12) 


In  this  computer  equation,  the  bendiiig 


moment  is  proport-ional  to  i(t)  d*'X  and  the  shear  is 


(it 


Z 


proportional  to  d 

dt 


i(t)  d^Z 
dt^ 


As  the  ship  acts  as  a  free- free  beaa,  the 
end  boundary  conditiorP  on  ths  prob2.eta  ar^   that 
the  bending  moment  and  shear  are  z^ro  at  each 
end ,  The  boundary  conditions  are  expressed  as 


d^I  ^  d 


i(o)  ilJ^  ==  - 
dt"^   dt 


i(o) 


d^X 


It 


«  i(T)  13  «  i 
,2 


dt^ 


-  0 


dt*"   dt  L 

The  computer  circuit  for  the  soltition  of  equation 
(12)  is  given  in  Fig«  16,  The  end  conditions . of  zero 
bending  moment  and  shear  are  sati?:^fied  on  the  rtcmputer 
by  initially  shorting  the  feedback  capacitors  of  k'^ 
and  A2  respectively.  As  the  slope  and  deflection  at 
each  end   of  the  beam  are  unknomij  these  are  simulated 
on  the  computer  by  battesy  voltages  -Va  and  Vb  x^espectively 
applicid  initially/  to  the  capacitors  of  A^  and  A50  ¥b 
was  fixed  at  about  6  Yoits  and  -Va  was  varied  for*  the 
different  trial  solutions  for  the  various  modes  until 
the  end  conditions  of  sero  shear  and  bending  moment 
were  satisfied. 

On  Table  I  is  tabulated  the  original  data  on  the 
APA  S79  the  ship  for  which  frequency  of  vibration  was 
desired «  This  ship  has  the  foliovjing  general  characteristics 
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L  «  Load  water  line  length  »  400' 

B  «  llolded  breadth  «  5^» 

D  »  Molded  depth  ==37' 

d  «  Full  load  draft  •«  15?  6" 

A  =  Full  load  displacement  -  6S00  tons 

E  «=  Young's  Modulixs  of  elasticity  '^  1,93  r-  XO- 

tonf/ft'^ 

Io«  Maximum  area  moment  of  inertia  '^  262^  ft^** 

J^'io'^  J-fextoum  total  mass  per  unit  L     ^-  2»0C6      g 

tons  .sec^ 

For  the  purposes  of  calculation,  th©  ship  was 
divided  into  20  parts  of  20  foot  lengths  each.     For 
each  section  there  is  tabulated  the  moment  of  inertia 
and  mass.     Fig,  17  sho?7s  the  distribution  of  mass  stA 
moment  of  inertia  of  the  APA  6?  as  introduced  into 
the  computer.     On  Table  II  there  is  tabulated  i(t) 
^.nd.  p{t)  %vhich  were  siraulated  by  iieans  of  40  steps ^ 
k  steps  ry^T  second,   on  the  stepping  relay  and  resistor 
panel,     ihe.  resistances  in  iiiegolims  added  for  each  two 
steps  are  tab-alatad  on  Table  IIo     The  negatibre .  sign 
before  the  :;'c:ji.iua:.ij;  -i  indicai^i^o   v.uai/  that 

•.■?e3is-Gance  wsis  subtract  >  f;:lve  t-he  p.  '  valua  of 

iCt)    o:;^  3(t)   for  that  step.     i-'Oi'  i(t)j  t/s  resistors 
listed  are  introd'aced  as  ijiput  resistance  to  A^^  ai^d  fv^r 
j^  ft),   the  resistors  listed- are  introduced  as  fsedoac^t 
resistance  to  Ai  to  satisfy  equatio-ii  (12) , 
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As  the  stepping  relay  resistor  panel,  Fig,  6, 
was  criginally  planned,  the  length  of  solution  on  the 
computer  should  have  been  T  •*  10  seconds.  However j 
as  the  stepping  relays  used  stopped  the  problem 
immediately  upon  reaching  step  40,  the  length  of  the 
solution,  T  *=  9*75  seconds.  This  condition  could  have 
been  corrected  by  rev/iring  the  stepping  relays*  control 
relay  circuit  to  provide  for  a  quarter  second  pause 
on  step  40  before  the  initial  condition  relays  again 
imposed  the  end  condition  and  stopped  the  problea* 
The  authors  of  this  paper  did  not  feel  at  liberty  to 
change  the  equipment  in  this  manner,  and  felt  that 
the  quarter  second  lost  could  be  accounted  for  in  the 
solution  knowing  that  T  was  actually  9»75  eeqonds. 
Even  thou^  the  quarter  second  in  the  length  of  the 
computer  solution  meant  ten  feet  in  the  ler,;gth  of  the 
ship,  this  "lost"  section  of  the  bov7  does  not  materially 
affect  the  vibratory  characteristics  of  the  vessel. 

Fig,  9  shows  schematically  the  arrangement  of 

the  computer  and  complete  network  of  controls  and  power 

supply.  Fig,  6  is  a  photograph  of  the  complete 

network.  The  outputs  of  A2  and  A^  are  shovv'n  connected 

through  amplifiers  to  the  tv^o  channels  of  the  Brush 

recorder  for  recording  oscillographs  of 

d^ 

ai^<i  -i(t)  £X  .  dt 

dt2     ' 


i(t)'d^X 
dt^. 


As  iD  thf,  solution  of  the  ur.ixoKi  free-free  besr<i, 

cor;.e«fc  solutions  she  ■•ulSlliaGrit  of  the 

end  conditions  of  zc^r  ■   - 

-uiiB  oi   aei  auii  oGndiog  moment  were 

obtained  vdien  the  minim  .^iniim  of  ^.i(t)   VK^ 

d^-omding  upon  the  nv.-' ■     'e  nod<P      -  ^  <-u 

e  i-ioae^    ,  I  through 

tne  aero  axis*     Th©     d    '~"•W•i-^    -i2sO 

in  measuring  the  lengta  o:f  the   solution,  T^ 

imh  the  eoinponents  of  the  computer  network 
oparating  as  descriDad  in  the  E  :ent  and  Pro 

sictio^i  of  this  paper,   correct   ......ions  w.x^  o^.aine. 

b;-  first  varying  ^7a  from  trial  to  trial  until  snd 
conditions  were  met.     Then  Q-ms  varied  from  trial 
'^o  trial^mtil  a  solution  of  length  T  -  9^75  seconds- was 
recorded.     It  was  found  that  with  each  trial  value  of 
0  inl^roduced,  «?a  had  to  be  changed  slightly  for  the  end 
coi^ditiors  to  be  satisfied  for  that  trial* 

^•he  oscillographs  of  solutions  obtained  fcr  the 
first  two  rAodes  are  given  in  Fig.   1^  and  Fig,  ISu 
it  VTas  ifi;'P03sible  due  to  the  limitationK  of  the  equip-^. 

^^^  *°  obtain  a  curr^cu  solu:tion  satisfying  the 

end  conditions  for  the  third  mode.     I---  r-^-^mpts  were 
2J£de  and  the  near  soluticns  obta'  otted  and 

cross  plotted  on  Fig.  20,     From  thin  fi^.^.^      :. 
approximate  value  of  G  n&oessary  for  a  correct  ,  i 

could  be- obtained.     This  method  of  plotting  varia  . 


:f?:7oa  correct  solution  waa  developed  by  Howe  in  Ref  ^ 
(6)  and  was  extensively  used  by  hiui  when  correct 
oscillograph  solutions  were  difficult  or  impossible 
to  obtain 4 

Theoretical  solutions  of  the  frequency  of 
normal  nodes  of  vibration  of  the  APA  ^7  have  been 
calculated  by  graphical  methods  and  by  a  calculator 
designed  by  IBM  ,  Comparative  results  are  listed  beloi*:, 


Frequency  of  Vibration  «•  Noiisal  Modes 
Binding    Only 


Mod© 

Computer 

IBM       1 

C 

rad/sec . 

rad/sec* 

1 

6*0 

12,4 

•£••*(>  a  P            1 

2 

0.9S 

29.^ 

2^.4      j 

3 

0,25 

60.6 

Considering  now  the  effect  of  shear  deflection 
and  rotary  inei^ia  in  addition  to  the  bending  de- 
flection effect  just  discussed: 

Equation  (7)  was  the  complete  differential 

equation  of  the  elastic  curve  in  vjhich  the  torm 

had  the  physical  interpr^^tation  as  a  mcment  du©  to 
elastic  deformation  in  rctary  inertia.  This  terra 
was  reduced  by  the  differential  operator  ^  "^ 


to  an  equivalent  distributed  load,  I /^  is  the  mass 
moment  of  inertia  of  a  cross  section  area  per  unit 
length  of  the  bfeara* 


In  equation  (7)  is  also  included  th^  tnTm.      EIa 

KAG 

which  is  also  reduced  by  the  operator  \3       to  an 

equivalent  distributed  load..  This  is  tb'r,   culy  tenst 
left  in  the  original  equation  (7)  representing  the 
effect  of  shear  deflection  as  the  tena 


-^-Za  i/  ^IXfe-^   ^^^  heen   left  out  in  settin 
"^t^  KAG  '^    "^  t^ 

up  the  computer  equation. 

Regrouping  the  teras  of  equation  (7)^  the  follow- 
ing fourth  order  differential  equation  is  obtained 
which  considers  the  "effect  of  bending  deflection  and 
an  approximation  to  the  effects  of  shear  de-flection  and 
rotary  inertia: 


:5 


l>'ii 


2 

2 


f 


«»>aKti'<«K-'*  r^    ^     -Ti-mta  ^^^    i       ■  ■-  ■  r^i^-iJirrfi  ~  V  -  t-i- 

'-  1)  x^  Ikag  a/  ^t^  '. 


(13) 

where  I  is  area  moment  of  5.nertia  of  a  cross  section 
area  ^ 

In  addition  to  the  assumption  that  yCx^t)  = 
X(3c)  ©'^''■^^s,  it  is  assumed  that  over  the  length  of  the 
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ship  x,ne   '^sriii/E^^  '*'  i  ^  ''"^  ^^  constan-fco  Sighting 
VKAG   A  / 

of  the  tabulated  data  for  the  ship  considered,  Table 
nijShovfs  the  reasonableness  of  this  assumptlona 
Denote  this  constant  as  B^  then 


'E 


Ikag 


A 


With  I  and  //   as  functicnn  only  of  x  dfstc 
the  beaia  and  independent  of  time,  we  have. 


d^ 


dx'^ 


^// 


a) 


2v 


0. 


long 


:i4) 


iS  was  done  for  bending  deflection  only^  l^t 

I  ^  I^i(^)  and^  ^//q^M   and  set  up  the  coDiputer 

equation  by  making  a   change  in  the  independent  "Variable ^ 

x^   in  the  original  equation  (14)  to  t,  ti),i@  in  seconds 

on  the  computer  solution.  Than  t   ^  T  sc^ 

L 

Making  the  necessary  substitutions  and  dividing  through 

&ya)^  andM^f  equation   (14)   bacomes 


111 


(l 


<a. 


Ci)  ^h^-  7(o  ^^^ 


Q 


i(t)   d% 


dt^. 


.J.  IqBT^  d^ 
l2  dt^ 


i(t)p(t)l   2 


3{t)I  =  0 


Lettins:  a 


ti 


(aJlm^^ 


G 


EI 


a 


,4 

2 


(15) 


J!? 


and  H  =  lo  B  T^  ,  equation  (15)  becomes 


l2 

G  d^  i(t)  d^X 
dt^  L    dt^  . 

^  H  d^ 
dt^ 

i(t)p(t) 

X  «  p(t)x  «  0  ■ 

The  natural  frequency  of  the  ¥ibration  in  the 

n"^^  mode  is 

- 

- 

,1/  Slq 

In  order  to  obviate  the  necessity  of  setting  up 
a  new  computer  circuit  to  introduce  the  product  term 
of  the  two  variables,  i(t)  and  ^(t),  it  vfas  felt  a  good 
approximation  could  be  attained  by  selecting  from  the 
data  available  a  representative  average  value  of  the 
product  term   i(t)  p(t)   which  could  be  assrued  con» 
stant  over  the  length  of  the  heem.     This  was  done  and 
X7ith  the  regrouping  of  constants  in  the  second  term  of 
equation  (16)  to  a  single  constant,  D, 

D  «  H  ri(t)  p(t)] 
and  the  computer  equation  becomes 


dt^L    dt^J     dt^ 


0. 


(17) 


The  end  conditions  for  a  free- free  bear-n  are  that 
bending  moment  and  shear  are  zero  at  both  ends.  In  the 
previous  problems  discussed,  these  were  proportional 
to  the  second  and  third  derivatives  of  the  deflection 
respectively.  IJhen  in  addition  to  bending  deflection. 


shsar  and  rotary  inertia  effects  are  considered,  the 


bending  moment  and  shear  force  are  given  as 


12 


M  B,  EI 


ur 


^^7^ 


KAG 


y 


Q  =  EI 


r 


1  ^/lL^    j^^.  /£id -f..^!^  J^I 


KAG 


'  Z)x^ 


L 


KAG 


E    'bx 


However,  in  setting  up  the  computer  equation 

the  assumption  iiras  made  that  _2)^  i^//.kIX 

dt^   KAG  ^t^ 

could  be  neglected  for  a  good  approximation  of  the 

effect  of  shear  and  rotary  inertia*  Following  this 

assumption,  the  bending  moment  is  proportional  to 


il 
V 


i(t)  S»^  and  the  shear  is  proportional  to 

.2 


d_ 

dt 


i{t)  ±4 
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in  the  computer  equation  (17). 

The  bD?indary  conditions  for  the  solution  of  the  computer 
equation  are  then  expressed  as 


iCo)^::! « i 


dt' 


dt 


[ 


i(o)  dfZ 
dt^J 


dt     dt 


i(T)^| 
2 


dt 


0 


The  computer  circuit  for  the  solution  of 

equation  (17)  is  given  in  Fig,  16,  where  the  dotted 

line  V7ith  resistance  C/D  from  the  output  of  A^^^  to  the 

input  of  Ag  is  included  to  accomplish  the  D  ^X  tona 

dt 


in  the  computer  equation.  The  procedure  for  the  so- 


lution  of  equation  (1?)  is  similar  to  that  described 
for  the  solution  of  equation  (12)  where  bending  de- 
flections only  were  considered «  In  addition  to 
varying  Gg  it  is  no^^j  necessary  to  vary  the  ratio 
O/l).     As  D  remains  constant  for  all  modes,  the  ratio 
1/D  eould  have  been  introduced  into  the  computer 
circuit  instead  of  C/D  and  ,5/C  would  have  been  the 
variable  introduced  as  a  feedback  resistor  on  A^j 
necessary  from  trial  to  trial  for  each  raode« 

The  CGcillographs  of  solutions  obtained  for  the 
first  four  modes  are  given  in  Figs.  21  through  24 
inclusive. 

It  was  found  that  the  addition  of  the  C/D  input 
resistance  to  A2  from  A/^  gave  considerable  stability 
to  the  computer  circuit  in  the  range  of  resistance 
used.  The  alternative  of  using  a  constant  1/D  ratio 
as  mentioned  above  as  an  input  resistance  vjas  not  as 
successful  in  stabilizing  the  circuit  and  so  was  not 
used , 

Theoretical  solutions  of  the  frequency  of  normal 
modes  of  vibration  of  the  APA  ^7  have  been  dete/r'mined 
by  a  calculator  designed  by  IW-l   for  bending  and  shear 
deflection  and  rotary  inertia  effects,   Goiaparative 
results  are  listed  below. 


ib 


Frequency  of  Vibration  -  Normal  Xlodes 

Bending  and  Shear  Deflection  and  Rotary  Inertia 


ilode 

Computer 

im-i 

1 

C 

rad/sec. 

rad/sec. 

7^00 

11.29 

10.30 

2 

1.6L^5 

23.02 

19.95 

■ 

3 

0«725 

35.00 

30.02 

4 



0.410 

46.60 

39. 2£^ 

The  frequencies  obtained  by  raeans  of  th@ 
analog  computer  v/hen  the  effects  of  bending  deflect iono 
only  are  considered  are  approximately  five  to  eight 
percent  higher  than  the  IM  solutions.  The  analog 
computer  results  are  calculated  directly  frora  the 
oscillographs  and  have  not  heon   corrected  fo:.'  recorder 
pen  lag  or  power  supply  frequency  variation.  Ths 
latter  is  an  important  factor^  correction  fo'f  which 
should  be  made  for  better  accuracy  of  results.  Line 
frequency  variation  has  an  effect  both  on  the  measure- 
ass  nt  of  -'T"  on  the  recorder  tape  and  on  the  balance 
of  ths  operational  amplifiers,  necessary  for  repetition 
of  solutions. 

Another  source  of  error  in  the  compute!^  effecting 
the  accuracy  or  the  results  lies  in  the  stepping  relays 
and  resistor  panels  used.  Every  effort  was  raade  to 
have  accurate  resistances  on  each  ster)  for  the  simu- 


lation of  mass  and  moment  and  inertia 


However 5  the 


J  J 


many  plug-in  connections  of  resistors  in  stacks  on 
the  resistor  panels  introduced  inaccuracies  in  the 
actual  resistance  obtained  for  each  step«  Then  too, 
it  is  known  that  the  bridging  contacts  of  the  stepping 
relays  did  not  always  perfectly  bridge  from  one  step 
to  the  next* 

It  is  felt  that  the  results  were  well  within 
the  accuracy  of  the  computer  network  itself,  and 
that  corrections  made  to  the  oscillograi:>h  records 
as  mentioned  above  would  improve  the  precision  of 
the  values  of  frequencies  of  vibration  obtained „ 

The  results  obtained  when,  in  addition  to  bending 
deflection  effect,  an  approximation  to  the  effects  of 
shear  deflection  and  rotary  inertia  was  considered ^ 
are  progressively  higher,^  (from  about  ten  percent  for 
the  first  mode  to  about  eighteen  percent  for  the  fourth 
mode)  than  those  obtained  by  the  IM  computer. 

It  is  apparent  that  these  deviations  which  increase 
with  the  higher  modes  result  from  something  more 
than  the  inaccuracies  in  the  computer  network. 
The  assumptions  made  to  approximate  the  effects  of 
shear  deflection  and  rotary  inertia  in  setting  up 
the  computer  equation  mv^st,  in  a  large  part,  account 
for  the  increasing  error «  These  assumptions  were 
necessary  to  avoid  complicating  the  present  computer 
network  to  a  degree  out  of  proportion  to  the  actual 
effect  of  shear  and  rotary  inertia  on  the  frequencies 
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of  vibration  of  the    ship* 

IJith  a  computer  constructed  of  more  precise 
and  stable  cc::iponents,   the  product  term  i  i(t)3(tjj 
could  be  introduced  as  a  variable,     A  netviork  could 
be  set  up  to  solve  a  computer  equation  which  includes 


^t2  KAG '     >)  t^ 

The  precision  of  the  values  of  frequencies  obtained 
from  such  an  analog  computer  netv/ork  should  be 
definitely  better  than  in  the  present  instance* 


CONCLUSIONS 


Solutions  to  many  engineering  problems  of 
practical  interest  involving  higher  order  dif- 
ferential equations  v/ith  variable  coefficients 
may  be  obtained  by  means  of  a  relatively  simple 
and  inexpensive  electronic  analog  computer^ 
Solutions  so  obtained  are  well  within  the  accuracy 
necessary  for  most  engineering  purposes. 

The  accuracy  of  solutions  obtained  are 
limited  by  the  precision  of  the  computer  com« 
ponents  used  and  regulatio2i  of  the  associated 
power  supplies.  The  assuznptions  made  in  reducing 
an  exact  differential  equation  to  a  computer  equation 
are  in  a  large  part  necessitated  by  the  precision 


'1 
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of  the  apparatus  used.  Low  precision  of  co.iponents^ 
for  instance,  M)uld  limit  the  numbei-  of  amplifiers 
and  variables  introduced  into  the  computer  network 
for  a  given  desired  accuracy «  As  the  coraplexity  of 
the  network  increases,  so  must  the  prec5.sion  and 
stability  of  the  components  inc3:'eas0o 

Instead  of  introdueing  another  variable  of  the 
original  equation  into  the  computer  ^fith  the  conse- 
quent additional  amplifiers  and  component  circuits, 
a  'Constant  average  effect  of  the  variable  msc  ""^ 
i':ntroduced.  Errors  resulting  from  such  asvS^^Mption 
of  average  effect  must  be  weighed  against  those 
resulting  from  lack  of  precision  of  the  circuit 
elements.  In  the  case  of  the  APA  B7   in  thie  paper, 
the  constant  average  effect  of  shear  deflection 
and  rotar]"  inertia  embodied  in  the  term  f  i(t)^(t)J 
vrou.ld  appear  not  to  be  representative  of  the  ti*ue 

effect.  However,  if  the  term  ^^     lu    jjYyiKit) 

>>t^   KAG    ^t2 

had  been  retained  in  setting  up  the  computer  equation, 

there  would  undoubtedly  have  been  close:.'*  areeemcnt 

in  results  obtained.  In  the  latter  case,  the  average 

effect  was  assumed  to  be  aero,  an  obviously  over 

simplification  in  light  of  the  results  obtained  by 

the  computer. 


Economy  of  time  and  personnel  is  the  primary 
benefit  of  such  an  ranalog  computer  as  described  in 
this  paper,  T/ith  the  necessary  equipment  available 
and  having  familiarity  with  the  operating  procedures, 
the  solution  of  higher  order  differential  equationsj 
vjith  variable,  coefficients  vjould  be  a  matter  of  a 
fmi  hours  for  a  single  operator. 

The  analog  computer  is  especially  adaptable  to 
the  solution  of  design  problesns  where  the  study  of  the 
effects  of  varying  design  parameter's  may  be  eondueted 
with  little  effort  by  simple  external  changes  to  a 
single  computer  netvrork. 
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CALCULATION  OF  FREQUENCY  OF  VIHaATIOK 
Bending  Onlv 


Constants  used: 


T  -  9.75  sec.      T^  =    95o0625  sec^ 


EI^         llT. 93x10^  X  262g  =^      0,3U 


'/V„L^  f   2.006  X  256  X  10^ 

1st  Mode 


C    =  6„0 
rV  -    "^^       -     95.0625     s    3S.a5 


yc  7^.0 


0),  -  o(, 


EIq       =:    3S.85  X  0,314    rad./sec, 


CO,  -  12,4  rad/sec, 
2nd  Mode 


C  =  Oo9S 


=  ^  =  95  >  0625  =  95ol 


X     _  cdf  ®^o        =   95a  X  0,314  rad./sec. 
<^t   =  29 a  8  r^dv/seco 


3rd  Mode 

C    =   Oo25 
o(  -  -1?    s     95.0625      s  190,125 

^  fcT    '   ]ro;25 


>i|  EI^       3.190.125x0.314    rad./sec. 
CO-i  ~    60.6  rsd./s©c. 
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CALCULATION  OF  FREQUSWCY  OF  VIBRATION 
Bending,    Shear,   and  Rotary  Inertia 

Constants  used; 

t"=    9.75  sec  T^   =    95.0625  see, 2 

SIq      =:   0.0986 

From  Table  III        B  =  1.30Z,    Vft^      j  [Cir)  /^(t)'^    r.  Q, 


45 


533 


^^   "     L^J^      "     2628  X  1„304  x  95  e  0625 
l2  160,000 

H    =      2c039  sec, 2 

D   =     H    [  i  (t)  B  (t)]     -    2.039  X  .533 

D   ^      loOa?  see. 2 


1st  Mode 


7oO 


a«    =   T^      =s   9040    =      1291 
i       — ==  " 


«>.   '  a^^    SIq     =    1291  X  0.09S6     ^  127.5 


<^j    =    11.29  rs.d./sec. 


2rid  jfod® 

C  =    1,685 

o:*  =■-   T-i-     «      9040    =         5360 


•£ 


C  1.685 


^e.^  ai  EIq    =  5360  X  0.0986    =  530 


^o^^ 


'^a  =•      23.02  rad. /sec, 
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3rd  KodQ 

C  ^=  Oo725 

a*  ^jt^'   29^   =  12,470. 
C    0,725 

OJ^    =  12470  X  0,0936   =  1232 

^3  =  35 oO  rad./sec« 


4th  ?4od® 
C      «    0.410 


as    =    t4   =,    9040      =    22,050 
'^      T'        0.41 


'^^^  ""        O'i   ?I2  "^    22050  X  0.09^6   «=    2175 
^-^  ==    46,6  racl./seco 


TABLE  I  APA  &7 


DATA  FOR  -CALCULATION  OF  NORMAL  MODES  OF  VERTICAL  VIBRATION 


400  ft. 


E  =  1.93  X  10°  tons/ ft 


Section 

Stem 
To  Bow 

I  ft> 

/4  ■  pt^ 

K 

^H? 

KAG 
tons 

a-i 

617 

0.1775 

0,310 

5o35 

lo27S 

l=-2 

1157 

0,6271 

0.271 

7ol$ 

lo494 

2-3 

15S6 

0.9025 

0.220 

9«30 

lc57S 

3-4 

1S95 

1.1300 

0,171 

12e,29 

I062I 

4-5 

2146 

1.2946 

Ocl50 

14.03 

1.622 

5-6 

2334 

1.3548 

0.149 

13«75 

lo579 

6-^7 

22;.  5  4 

1.456g 

0J.5(> 

13.61 

1.637 

7-8 

2532 

1.9493 

0.161 

13.82 

1«715 

g-9 

25^^5 

1.9540 

0.156 

14.44 

1.737 

9-10 

2609 

2 0 0060 

•   O.li/ 

14.93 

1,658 

10^11 

2623 

1.9401 

0.137 

14.65 

1.548 

ll-^l?. 

2626 

1.9159 

0,134 

14o93 

1.542 

12-13 

262s 

1.6765 

0.141 

15.97 

lc736 

13»14 

2614 

1.4242 

0.152 

15.49 

lc815 

14-15 

2599 

1.2439 

0,16S 

12.99 

1.682 

15-16 

2A.93 

1,0139  . 

0,186 

12.50 

1.792 

16-.17 

22S1 

1  ■ 
0.7112     !     0.205 

llo30 

1.S66 

17-13 

1934 

0.448S 

0.225 

10.14 

1.759 

18-^19 

1447 

0.254c^ 

0.245 

7c  99 

1.309 

19«20 

738 

0.2613 

0,264 

3o49 

I0II7 

TABLE  II     APA  87 
i(b)   and/S(t)    in  Megohms  as  Introduced  into  Computer, 


I^  «  26^8  ft 
o 


yU  Q  «  2  0060  tons  sec'/f't'* 


X  =  lo  i(t) 


Stern  to 
Bow  Section     i(t) 


AR  MEG  AR  J'^EG  ^^  ^^  ^^  ^^  ^^ 

Ai(t)  Breakdown  /Bit)         A/6(t)        Breakdovm 


DATA  FOR  CALGULATIOH  OF  CONSTAHTS  B,  H,  and  D . 

^kIg  -^  i  )  »   H  =  lq£r ,   D  --=  H  [i(t)  /5  (t)] 


B 


Section 

Vkag 

A 

KAG"*"  A 

i(t)/3(t) 

0-:' 

lo511 

0.187 

1.698 

0.022 

l'-2 

1.291 

0.140 

1.431 

0,137 

2«3 

1.222 

0.107 

1.329 

0.272 

3-4 

1.190 

0.081 

1.271 

0.405 

4-5 

lolgg 

^    0.071 

1.259 

0.,526 

5-6 

1„222 

0.073 

1.295 

o„6co 

6-7 

ij.eo 

0.073 

1.253 

00677 

7-8 

1.124 

0„072 

1„196 

0.935 

g->9 

1*111 

0.069 

laso 

0.959 

9-10 

lcl63 

0.067   ' 

1.230 

0.993 

10-11 

1,247 

0.063 

lo315 

Oo965 

11-12 

1.250 

0.067 

1.317 

0.955 

12-.13 

1.111 

O0O62 

lcl73 

0,836 

13-14 

lc062 

0.065 

1.12'^ 

Oo706 

14-15 

1.U7 

0.077 

1,224 

0,615 

15-16 

lc076 

0,030 

I0I56 

0.481 

16~17 

lo033 

0.085 

1,118 

0.3C8 

17-1.1 

lc097 

0.09s 

1.195 

0<,171 

18=.19 

1,279 

0„125 

lo404 

0.073 

19-20 

1.730 

0.1S2 

lc912 

C„036 

£=  26.083                   t  =  10.671 

5  =     1-304              ^^E-°     0.533 
Ave.. 
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Figure   7 
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Figure  11 
Oscillograph  Solutions 
of  Deflection  of 

Clamped  End  Beam. 
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Figure  13 
Oscillograph  Solution  of  Deflection  of  Hinged  3nd  Beam. 
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Figure  15 
Oscillograph  Solutions  of  First  Three  Normal  Modes, 

Uniform  Free -free  Beam. 
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